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Abstract 



Using new scaling parameters (3i, we derive simple expressions for the ex- 
cess thermodynamic properties of the Mean Spherical Approximation (MSA) 
for the ion-dipole mixture. For the MSA and its extensions we have shown 
that the thermodynamic excess functions are a function of a reduced set of 
scaling matrices F^. We show now that for factorizable interactions like the 
hard ion-dipole mixture there is a further reduction to a diagonal matrices 
(3^. The excess thermodynamic properties are simple functions of these new 
parameters. For the entropy we get 

(kV_ 



5 = -|^|(^[/3.])„ex 

where is an algebraic functional of the scaling matrices of irreducible rep- 
resentations X of the closure of the Ornstein-Zernike. 

The new scaling parameters Pi, are also simply related to the chemical 
potentials of the components. The analysis also provides a new definition of 
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the Born solvation energy for arbitrary concentrations of electrolytes. 



Typeset using REVT^ 



ACS 61.20.Gy. 



I. INTRODUCTION 



The remarkable simplicity of the mean spherical approximation (MSA) [Q-^ and its ex- 
tensions using Yukawa closures can be summarized by the fact that for a wide class 
of systems the entropy has the same simple functional form. The MSA |^ is the solution 
of the linearized Poisson-Boltzmann equation, just as the Debye-Hiickel (DH) theory 0. It 
shares with the DH theory the remarkable simplicity of a one parameter description [p!0| ( 
the screening length n) of all the thermodynamic and structural properties of rather diverse 
systems. The major difference is that in the MSA the excluded volume of all the ions are 
treated exactly. 



The MSA thermodynamic properties of real electrolyte solutions are expressed by sim- 
ple analytical formulas which are in many cases remarkably accurate. This is so because 
the MSA satisfies a number of exact asymptotic relations, such as the large charge, large 



density limits of Onsager ||TlHl4| and the large charge, small density limits implied in the 
Wertheim-Ornstein-Zernike equation. Very recently |^,^we have been able to extend the 
MSA closure analytical solutions to any arbitrary closure that can be expanded in damped 
exponentials (Yukawa functions), and obtain explicit, analytical forms of the excess thermo- 
dynamic functions in terms of a matrix of scaling parameters ( the EMSSAP, or Equivalent 
Mean Spherical Scaling Approach [0). 



For systems with Coulomb and screened Coulomb interactions in a variety of mean 
spherical approximations (MSA) it is known that the solution of the Ornstein Zernike (OZ) 
equation is given in terms of a single screening parameter F. This includes the 'primitive' 
model of electrolytes, in which the solvent is a continuum dielectric, but also models in 
which the solvent is a dipolar hard sphere, and much more recently the YUKAGUA model 
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of water [|T^] that has the correct tetrahedral structure. The MSA can be deduced from a 



variational principle in which the energy is obtained from simple electrostatic considerations 
and the entropy is a universal function. For the primitive model it is 

AS = -kV— 

where F is MSA the screening parameter. In general it will be of the form AS = 5'(r) which 
is independent of the geometry of the problem. F is a general, non diagonal, scaling matrix. 
We have shown that in all known cases the scaling matrix F is obtained from the variational 
principle 

Ionic solutions are mixtures of charged particles, the ions, and the neutral solvent par- 
ticles, most commonly water, which has an asymmetric charge distribution, a large electric 
dipole and higher electric moments. Because of the special nature of these forces the charge 
distribution around a given ion and the thermodynamics does satisfy a series of conditions 
or sum rules. One remarkable property of mixtures of classical charged particles is that 
because of the very long range of the electrostatic forces, they must create a neutralizing 
atmosphere of counterions, which shields perfectly any charge or fixed charge distribution. 
Otherwise the partition function, and therefore all the thermodynamic functions, will be 
divergent [|1^]. The size of the region where this charge shielding occurs depends not only 
on the electrostatics, but also on all the other interactions of the system. For spherical ions 
this means: 



1. The internal energy E of the ions is always the sum of the energies of capacitors. For 
spherical ions the capacitor is a spherical capacitor, and the exact form of the energy 
is 



where z* is the effective charge, (3 = 1/ kT is the usual Boltzmann thermal factor, e is 
the dielectric constant, e is the elementary charge, and ions i have charge, diameter 
and density ezi, ai, pi, respectively. For the continuum dielectric primitive model 
Fj = r for all i. 

2. The Onsagerian limits. When the ionic concentration goes to infinity and at the same 
time the charge diverges, then the limiting energy is bounded by 

AE = --J2p,zX (1.3) 

obtained by setting Fj ^ oo 

3. A further exact limit is the DH limiting law, which simply requires that for all ions in 
the system 

2Ti^K with = ^I^f2p^zl (1.4) 

4. Finally in systems that are strongly associating in the limit of total association the 
above equation still holds. This means that if component 1 forms a n-mer the DH 
limiting law must satisfy 

47r/?e2 



m 

i=2 



1.5) 



This limiting law is satisfied to lowest order in the multipole expansion by the closures 
of the Werteim-Ornstein-Zernike equation [p!7| , p!8[ . 

The ion-dipole mixture is the simplest model of a discrete solvent electrolyte which has 
analytical solutions |jl9|]- [^. The solution and thermodynamics are given in terms of three 



parameters bQ,bi,b2, which correspond to ion-ion, ion-dipole and dipole-dipole interactions, 
or their corresponding scaling parameters F, A, B, which are functions of only two coupling 
parameters for the ions and the dipoles do = UDebye and d2 = Si/dipoie (see Eqs.(|1.7|),( |1.12| )). 

In this paper we show that introducing two new scaling parameters (3o,(3i the problem 
is reduced to only two equations, which are easily solved by iteration. Furthermore, we get 



a rather compact expression for the excess entropy. The chemical potential of the ions and 
the dipoles is directly related to the scaling parameters, as we will see below. 

This also leads to a natural extension of the Born cycle to arbitrary ionic concentrations, 
and hence a robust general definition of the Born solvation energy for arbitrary concentra- 
tions. 



The MSA and its extensions have been used in many very interesting applications ||35| , |36 



Especially noteworthy are the engineering applications of S.Lvov and his group, and the 
Beijing group |]39|,|38[|. 



In section 1 we summarize previous work on the electrostatic equivalent models for ions 
and for dipoles. In section 2 we review formalism for the ion-dipole mixture using the old 
scaling parameters F, A, B. In section 3 we use a matrix formalism to introduce the new 
scaling parameters. In section 4 we derive the new results for the thermodynamics and the 
Born solvation thermodynamics. 



A. The electrostatic equivalent models 

The MSA can be interpreted using an electrostatic analog from which all the results of 
the MSA can be deriveded ||I^j40|j4T[| . For charges this is simply a spherical capacitor. Here 
we include also the dipolar case, for which the equivalent capacitor is a dielectric sphere 
surrounded by a media of different dielectric constant. 



1. Charge- Charge Interactions 



For the primitive model of ionic solutions in the general case [^the parameter F is 
determined from the equation 



where tho ionic charge is ZiC and number density pi = N'i/V, where A/) is the number of 
ions and V is the volume of the system. Also k, is defined by Eg. ( |1.5| ). We remember that 
the excess internal energy is 



[I. 



which is simply the energy density of a collection of spherical capacitors. 
The excess entropy is 

OTT 

Then F is determined in every case by the simple variational equation ( |1.1| ) 

d\j3AE{T) + rV(37r) 



dV 



(1.10) 



The functional form of the entropy is dictated by the excluded volume property of the MSA. 
Not every functional will guarantee no overlap of the hard cores of the ion pairs. 

2. Dipole-Dipole Interactions 

dddd 

For a system of hard spheres with a permanent dipole moment fi the MSA result can be 



expressed in terms of a single parameter A . Following Wertheim we have 



where 



4= \ M l - — ) (1-11) 

9 ew 
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^2 = ^ = (1-12) 



and ps is the solvent number density, y is the adimensional Debye coupling parameter. The 
MSA dielectric constant ew is given by 



..^^ (1.13) 

As has been often done in the hterature, the parameter A can be computed directly from 
the dielectric constant ew using the above cubic equation. This parametrization defines 
an effective polarization parameter. Just as in the case of the ions, there is a physically 
meaningful way of interpreting the MSA for point dipoles using the variational principle 
( |1.1| ). The dipolar system can be represented by a collection spheres with a dipolar charge 
distribution |T^: The parameter A can be interpreted in terms of a model consisting of 
a point dipole immersed in a sphere of dielectric constante^n which is surrounded by a 
continuum of dielectric constant eout- We write 



A = 

^out 

Now the MSA solution is given by the excess energy parameter 62 

60 A - 1 



e// 
9k 



6 A + 2 

where 62 is the dipole-dipole energy parameter defined below Eq. (|2.6| ). Since this expression 
is the Clausius-Mossotti equation we can interpret g^J^ as the effective Kirkwood parameter 
for this system. From here we calculate the induced dipole 

since d2 is the dipole of our spheres in reduced units. The excess energy 

(3E = d2Xd (1.14) 



The closure equation ( |1 . 1 1| ) can be rewritten as 



(A + 2)2 (A + 1)4 

This corresponds to exactly Eq.( |l.l| ) in the form 



A^-TV^ (1-15) 



d^E _ n dS _ , 16 

Vd\ VkdX (A + 1)4 ^ ' ' 



which can be integrated to yield 



Vk 3 



A + 1 



- 1 



[1.17) 



Now if we define the scahng lengths for the irrep x = 



To = A 



and for % = ±1 



A + 1 



then 



TT 



-S = - 

Vk 3 



- 1 



(1.18) 



Notice that they satisfy the Wertheim 'density' of the irreps since they are obtained by 
setting 



pi = -(l/2)po 



Furthermore, the choice 



^out 1 ^in A 



which corresponds to placing the dielectric sphere with an embedded dipole in the vacuum 
satisfies all the physical requirements of the case. 

II. lON-DIPOLE INTERACTIONS 



.1. In this section we will deduce new equations for the MSA thermodynamics for 
ion-dipole mixtures. 
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A. Summary of previous results 



For the sake of clarity we summarize former work [p!9|,p0|,p^,p^,p5|-p7| . Expanding the 



total pair correlation h{12) in rotational invariants [53] 



h{12) = /i°°°(ri2) + /i""(ri2)$"" + h'''\n2)^''' + h'''>{n2W + h''''{ru)<l>'''' (2.1) 



^011/^ \^on I £,101 



101 I £110/ 



,110 I £,112/ 



where /i'""^(ri2) are the coefficient of the invariant expansion, which depend only on the 
distance ri2 between spheres 1 and 2. The rotational invariants <|)™-"^ depend only on the 
mutual orientations of the molecules. The relevant correlation functions are 

• ion-ion: 



hiiir) 



(1/2) 



000, 



r) -/i1°°(r) 



(2.2) 



MO = (1/2) hl'}{r)-h'_'}{r) (f-A) 



(2.3) 



dipole-dipole: 



hdd{r) = -V3hYdir)fii ■ fi2 



+ Yy^dd^('^) [3(f ■/ii)(f ■/i2) - Ai -h] 



(2.4) 



where fi is the unit vector in the direction of /i. The solution of the MSA is given in terms 
of the 'energy' parameters 



lon-ion: 



ho = 27rp, / drhT{r)r 
Jo 



(2.5) 
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ion-dipole: 



b, = 2nJP^fdrhJ{r) 



(2.6) 



dipole-dipole: 



2 /-oo h^^^(r) 



(2.7) 



which, as will be shown below are proportional to the ion-ion, ion-dipole and dipole-dipole 
excess internal energy . These parameters are required to satisfy the following equations 

m 



2 I 2 j'A 



(2.8) 



aiKio - a2[l - Kn] = (ioc?2 



(2.9) 



(2.10) 



where do and d2 are defined by Eqs.(|1.7|)and ( |1.12| ). Also 



V = l + B 



(2.11) 



with 



B 



hi 6?(A + 2) 
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(2.12) 



/3a = l-| 
6 

In our previous work simple sets of equations were obtained using the scaling lengths p5|j2( 
r, A and B. These are related to the excess energy parameters and 62 of Eqs. (|2.5| - |2^ ): 



-r 

i + r 



+ B 



1 1 

+ 



i+r 2+ A 



(2.13) 
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2 + A 



(2.1 



and 



62 = 6 



A-1 
2 + A 



(2.1 



Then the closure equations |2.8|-|2.1CI| are: 



Ion-ion 



do 



{do^eJ{2 + 2V)-V} 

(i + A)(r + A) 



(2.1 



Dipole-dipole 



{yl + dl) _ 1 (^,^ BT 



(2 + A)2 9P V (2 + A)^ 



(2.1 



with 



= (1 + r + A)' (3 + r + A)' - ^(1 + r)'(i + A) (5 + 2r + 3A) (2.1 



and 

lon-dipole: This equation is in its original form 

(T+T) ^ 



(i^T^xf (3 + r + A)-(i±I)(i±4ii±lI±M 



(2.1 



which is really equivalent to the linear explicit equation for B 



B 



1 + r + A 



-sr (i + r) + doA(i + A)' 
rfo A (1 + A)' + 8 (1 + r) A (1 + r + A). 



(2.2 



so that our problem is reduced to a set of two equations for the two scaling paramet 

r,A. 
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Using Eq( |2.16|) we write 



with 



y 



1 - 



V (l + T + X) 

s(i + r)(i + A) V 



y 



(2.21) 



vrx 

For infinite dilution we get the hmit 

B 



1 



A (1 + r + A) 



K - 1] 



do 



2(i+r) 



f2.22) 



+ 

The interesting remark is that using the definition of B, Eq.( p.l2D , we get for the ion-dipole 
excess energy parameter bi Eq.(p.6|) 

3 ' 



bi = do- 



1 



1 



(2.23) 



'(1 + A)(2 + A) V V e«,. 
An efficient iterative solution of the equations is to start assuming that B = 0, and compute 
A 



1. Then solve Eq. ( |2.17| ) by iteration. Since A > 1 



A = -1 + Wl + 



3(^2 



(2.24) 



with 



C 



\ 



1 1 B 



-1 + 



A2 (2 + A)^ 



(2.25) 



where T is defined by Eq. ( p.l8|) . 



2. We get the physical branch of the ion-ion interaction equation (|2.16| ) 

-1 + 



2 



1 + 2- 



1 + 



-2r(i + r) + ^oTi;:; 
2(i + r)(i + A)(r + A) 



(2.26) 



(2.27) 



We illustrate the results of this procedure with a numerical example in which the solvent is 
'dipolar' water of dielectric constant = 78.4, density p = 0.0332 and diameter a^j = 2.7QA. 
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FIGURES 



FIG. 1. The screening MSA parameter F for the ion-dipole mixture compared to the primitive 
model Fo which assumes a dielectric constant eyj = 78.4 



FIG. 2. The polarization MSA parameter A for the ion-dipole mixture as a function of the 
molar concentration c 

The dipole moment of the 'water' is fi — 2.21 Debyes, which is in ther range of pretty 
much any other model of water. For the time being we discuss only the restricted primitive 
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model where everybody has the same diameter, which ideally would correspond to a solution 
of NaF. The general case of arbitrary sizes is left to future pubhcations. 



FIG. 3. The dielectric constant €w as a function of the molar concentration c 



FIG. 4. The reduced ion-dipole MSA parameter H as a function of the molar concentration c 
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In Fig. 1 to 4 we show the results for the MSA parameters F, A, B, and the dielectric 
constant as a function of molar concentration of a 1-1 salt for the restricted primitive 
model in which the diameters are all equal to 1. 

1. The scaled relations 



Using the scaling results [3,P^,p5|-p7| we find the new set of relations 



Also 



Q 



id 



2Df 



V 



V 



Qdd — 2 3 



1 + bo _ 
2Df \V 



{-A + S(l + F + A)} 



/2 ^ 



(2.28) 
(2.29) 



a.= (^){i + r}{F + A} 



(2.30) 



Q'm=[^{1-X' + B{2 + T + \)} 



(2.31) 



with 



2(1 + Fa) 



'hi2 + xy 

6 



2(1 + Fa) 



(2.32) 



B 



'h{2 + \y 

6 



(2.33) 



g' = 2 - 8- 



(2.34) 
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III. MATRIX REFORMULATION 



To obtain further simplifications we need to rewrite the solution of the ion dipole MSA 
1^,22] in matrix form. Remember that 



Q(r) = (r - 1)Q' + (l/2)(r - 1)^Q" - A 



(3.1) 



so that 



Q(0) = -Q' + (1/2)Q" - A; Q'(l) = Q'; Q'(0) = Q' - Q" 



(3.2) 



with 



Q' 



Q'ii Q'id 

Q'di Q'dd 



(3.3) 



Q" 





Q" di Q" dd 



(3.4) 



A 



ai 02 




The factor function matrix at the origin is 



Q(o) 



We define first the matrix 



-Q'ii - cti -Q'id - '^2 

-Q'di + Q" dif^ —Q'dd + Q^ddf^ 



with the inverse 



2D, 



1 + &0 &i/6 
&i/6 /^e/S 



/56/3 -6i/6 
-6i/6 l + 6o 



(3.5) 



(3.6) 



(3.7) 
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We get 



where 



and 



A = -2 



1 




1 
2 



Q(0) = 21 - T>f = 2Q° + 



where I is the 2 dimensional unit matrix and 



We use the well known relation 









^ id 

2 






Q^dd 




2 





1 _ 

2Dj 




2 



bi 9 _ 3(l+bo) 
4Df ^ 2Df 



(3.8) 



(3.9) 



(3.10) 



Q'(0) + [Q'(0)]^ = -Q(0) ■ Q(0) 



(3.11) 



to find the matrix Q'(l) (we use a — 1 throughout) which is also the matrix of the contact 
pair correlation functions: 



Q'(i) 



-2 
2 



no 9!d 
go 



no Qli 

^ii 2 

Oil 

2 ^dd 



2 ^dd 

This immediately suggests the new scaling parameters /?o, /3i 




2 



§0 
2 



with 



(3.12) 



(3.13) 



/3o = (1 + bo)do - hid2 
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A = — ^ + 



Then, we get 



Q'(i) 



-1/2(31 PoPi 



-PoPi -2 + 2/52 
From Eq. (|3.12|) we see that 



+ yt 



-61V2 



-blyl + ^[bl/A+{bo(3,-bl/12r] 



1/2 



and 



Pi 



^v+^'yv (l+A)^ 



-Ph! + ^Ibl/lQ + ((1 + bo)P3/2 + bj/Uf] 



1/2 



\ (1 + A) 
Furthermore, we get the new cross relation 



2 _ 16 B{1 + F) (1 + r + 2 A) 



V 



[b^Po/A-P,P,f = bl/A + Pl~yl^^ 
We observe also that the closure relations Eqs( p.8[ - p.lO| ) are equivalent to 



2,2 t2 

0-1 + Ci2 = % 



aib 



- - 02/53 = do A A 



1 2 



b_i 
2 



where we have used the definition 



+ Pl = A'A'' + ylA' 



A = do- 



Mo 

4A 



I + &0 + 



P& 



Pe 
A 



^Vb-p, 
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IV. THERMODYNAMICS 



We start with the equations derived by Vericat and Blum : The internal energy is 



(3E/V = — hodo^ - 2bido 62 - 2h2d\ 

4:71 L 



(4.1) 



which can be written as 



f3E/V = — [doiPo - do) + I2d2{(3i - d2)] 

4:71 



(4.2) 



For the Helmholtz free energy we got 
1 



PA/V = (-2 bo do' + 2hdo d2 + 2q'' + Q'J + 2 Q',/ + [Q[^f) (4.3) 



FIG. 5. The excess internal energy [3{E{c) — E{{S))/V as a function of the molar concentration 



From here the entropy can be computed using 

S (3 



which leads to 



^ = (^ot^o' - 4 hM - Qb2dl + 2q" + Q'J + 2 Q'J + (4.5) 
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Using the results of last section we get 



kV 
^ " 3^ 



16 V 32 



(1+Ar (i+A) 



3 ~ ^ 



16(A-1) fl IQV 



[i + xy \2 (l + Ar 



(4.6) 



where 



(4.7) 



FIG. 6. The excess entropy (3{S{c) — 5(0))/!^ as a function of the molar concentration c. 



The ionic excess energy as a function of the molar concentration is shown in Fig. 5. The 
excess entropy is shown in Fig. 6. 



This result can also be expressed in terms of the regular scaling parameters F, A, B [ p5 
16 P 



with 



3 + 



32 



;i + A)^ (i + A)= 



-F. 



1 16 (-1 + A) (I 



16P 



2 (1 + A)^ 



+ F 



(4. 



p I p 16P , 2ig(l + F) (1 + A) 



;i + A)^ 



V 



(4.9) 
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A. The activity and osmotic coefficients: An extended Born solvation formula 



Using the expression of the Helmholtz free energy Eq. (|4.3| ) we derive the chemical po- 
tentials of the ions 



(4.10) 



and of the dipoles 
Then the Gibbs free energy is 



{d(3A/V) 






_ (3e^ ' 




dpi 




do 




do 


{df3A/V) 






A(3p' 




dpd 


4:T!-Pd 


."2 


a3 


d2 



pG 



j3E 

— = pi(3pi + pd(3pd = 



(4.11) 



(4.12) 



where we used Eq.(^ 



The new scaling parameters 'diagonalize' the internal energy, since the internal energy 
consists now of only two terms 



E/V 



do 



- 1 



Pi 

d2 



1 



1 

V 



(4.13) 



The first term clearly is the charge energy term and the second term is the dipolar 
contribution, which however is dependent on the ion concentration. Since they are strictly 
additive one can construct a Born cycle for any arbitrary concentration. 

From Eq. (|3.15|) we get a general the Born solvation energy ||45|1 . 

^A^/£(l + 3^) 



E 



Born 



g(i+r)(i+A) 

V 



(4.14) 



[(1 + r) (fa 

where 3^ is defined by Eq.( p.22D . The dilute limit expressions for &i, the ion-dipole excess 
energy parameter were given by Blum, Chan et al and Zhou et al In all of these 

the solvent polarization contribution discussed by Kusalik and Patey is missing. Our 
current definition includes this contribution. Taking the zero concentration limit of Eq. ( ^T4D 
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we get the limiting Born solvation energy for a molecular solvent, which agrees with previous 
results as can be seen from equations ( |1.11| ) and ( [2^ ). We get 



/3/x,(0) 



560.0 



An 



1 



(4.15) 



.l + Ao_ 

where the subindex indicates the pure solvent parameters. A simple approximation is 
obtained when the concentration dependent A and are used: 



n papprox 
Born 



560.0 

(TiA , 



1 + A(a 



1 - 



(4.16) 
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We compare in Fig. 7 the exact expression and the approximate low density one, but 
with the correct density dependent A and e^. 



FIG. 7. Comparison of the MSA exact result for the Born solvation energy Eg. ( [4. 14 ) to the A 
approximation Eg. ( [4. 16 ) 



The activity coefficient of the salt referred to infinite dilution is 
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FIG. 8. The ionic activity coefficient log^± as a function of the molar concentration 



ln7± 



560.0 



rA./£ (1 + 3^) 



1 



l + Ao 



A' 



1 + r) (r A - ^ii±§(i±^ 

Fig. 8 shows the ionic activity coefficient as a function of the molar concentration. 
Similarly, we get the corresponding expressions for the solvent using Eg. ( p.l6D 



(4.17) 



with 



3c| 

T^Pd 



2 + A 



:i + r) (i + r + 2A) + p (i-^ 



1 + 



_9BW 

A2 (2+A)2 



V 



n 



:i + r + A)'(3 + r + A)' 



:i + rf (i + Af (3 + ^ 



9P 



(4.18) 



FIG. 9. Comparison of the MSA exact result for dipolar solvation energy Eg. ( 4.1^ ) to the A 
approximation . ) 



and the zero concentration limit is 



/3/i.(0) = ^ 



2 ^ 3 



2 + An 



.%er g I mol 



2 + An 



T = 298.13°ir 



(4.19) 
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which corresponds to the polarization energy Eg. (|1 . 14|) discussed in section ??. 

A simple approximation for highly polar solvents like water is to use the full concentration 
dependent A in these expressions. Comparisons are shown in Fig. 9 and 10 for the ion 
hydration and solvent hydration fromulas. 



FIG. 10. The dipolar activity coefficient logjd as a function of the molar concentration 



We remark that in the MSA the relation, 



G = E 



(4.20) 



still holds, and our chemical potentials satisfy the Gibbs-Duhem relation, 
The excess pressure can also be computed The expression is |^| 



P/ksT = S/VksT 



(4.21) 



The extensions to the general arbitrary size mixtures, as well as applications will be 
discussed in a forthcoming publication. 



.2. 
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